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Spontaneous symmetry breaking occurs in a system when its Hamiltonian possesses a certain symmetry, 
whereas the ground state wave functions do not preserve it. This provides such a scenario that a bifurcation, 
which breaks the symmetry, occurs when some control parameter crosses its critical value. It is unveiled that 
the ground state fidelity per lattice site exhibits such a bifurcation for quantum lattice systems undergoing 
quantum phase transitions. The significance of this result lies in the fact that the ground state fidelity per lattice 
site is universal, in the sense that it is model-independent, in contrast to (model-dependent) order parameters. 
This fundamental quantity may be computed by exploiting the newly-developed tensor network algorithms on 
infinite-size lattices. We illustrate the scheme in terms of the quantum Ising model in a transverse magnetic field 
and the spin 1/2 XYX model in an external magnetic field on an infinite-size lattice in one spatial dimension. 
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Introduction. Quantum phase transitions (QPTs) 11112] arise 
from the cooperative behaviors in quantum many-body sys- 
tems, in which long-range orders emerge. In the conven- 
tional Landau-Ginzburg-Wilson paradigm, the most funda- 
mental notion is spontaneous symmetry breaking (SSB), with 
the symmetry -broken phase characterized by the nonzero val- 
ues of a local order parameter. An SSB occurs in a system 
when its Hamiltonian enjoys a certain symmetry, whereas the 
ground state wave functions do not preserve it fsl Ql- The 
implication of an SSB is two-fold: first, a system has stable 
and degenerate ground states, each of which breaks the sym- 
metry of the system; second, the symmetry breakdown results 
from random perturbations. This leads to such a scenario that 
a bifurcation, which breaks the symmetry, occurs, when some 
control parameter crosses its critical value. Conventionally, 
this is reflected in local order parameters. 

The latest advances in our understanding of QPTs origi- 
nate from theperspectives of both entanglement [5] and fi- 
delity iS B B SIj which are basic notions in quantum in- 
formation science. In Refs. it has been argued that 
the ground state fidelity per lattice site is fundamental in the 
sense that it may be used to characterize QPTs, regardless 
of what type of internal order is present in quantum many- 
body states. The argument is solely based on the basic Postu- 
late of Quantum Mechanics on quantum measurements, which 
implies that two non-orthogonal quantum states are not reli- 
ably distinguishable 1.10,1 . In other words, the ground state 
fidelity per lattice site is able to describe QPTs arising from 
an SSB and/or topological order 111 ill . This has been fur- 
ther confirmed in Refs. 1 12, 13t], where topologically ordered 
states in the Kitaev model on the honeycomb lattice and the 
Kosterlitz-Thouless phase transition are investigated from the 
fidelity perspective, respectively. Moreover, even for systems 
with symmetry-breaking orders, it is advantageous to adopt 
the ground state fidelity per lattice site instead of using the 
conventional local order parameters, due to the fact that it is 
model-independent, although one may systematically derive 
local order parameters from tensor network (TN) representa- 



tions of quantum many-body ground state wave functions by 
investigating the reduced density matrices for local areas on an 
infinite-size lattice il4ll . However, it remains unclear whether 
or not it is possible for the ground state fidelity per lattice site 
to capture bifurcations arising from an SSB. 

In this Letter, we attempt to fill in this gap. First, we demon- 
strate that the newly-developed TN algorithms on infinite-size 
lattices may produce degenerate ground states arising from 
an SSB, each of which results from a randomly chosen ini- 
tial state subject to an imaginary time evolution. Second, it 
is unveiled that an SSB is reflected as a bifurcation in the 
ground state fidelity per lattice site for quantum lattice sys - 
tems undergoing QPTs with symmetry -breaking orders lll5ll . 
The significance of this conclusion lies in the fact that, on 
the one hand, this establishes the connection between the fi- 
delity approach to QPTs and the singularity theory; on the 
other hand, it is of practical importance since it makes pos- 
sible to locate transition points without the need to compute 
the derivatives of the ground state fidelity per lattice site with 
respect to the control parameter In contrast, the von 

Neumann entropy, a bipartite entanglement measure, fails to 
distinguish degenerate symmetry-breaking ground states. We 
illustrate the general scheme in terms of the quantum Ising 
model in a transverse magnetic field and the spin 1/2 XYX 
model in an external magnetic field. Here, it is worth em- 
phasizing that, although the scheme is applicable to quantum 
lattice models in any spatial dimensions, we restrict ourselves 
to quantum systems on an infinite-size lattice in one spatial 
dimension. This is achieved by exploiting the infinite matrix 
product state (IMPS) algorithm initiated by Vidal |l7l]. The 
extension to quantum lattice systems in two and higher spa- 
tial dimensions, which requires to use the infinite projected 
entangled-pair state (iPEPS) algorithm OlSll . is deferred to an- 
other publication [19]. 

Infinite matrix product state algorithm and spontaneous 
symmetry breaking. For quantum many-body systems on an 
infinite-size lattice in one spatial dimension, Vidal 1 17] has de- 
veloped a variational algorithm to compute their ground state 
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wave functions based on their MPS representations, which is 



a variant of the MPS algorithm Ii20l 12111 on a finite-size lat- 
tice in one spatial dimension. Here, we briefly recall the key 
ingredients of the algorithm. Assume that the Hamiltonian is 
translationally invariant, and consists of the nearest-neighbor 
interactions: H - 2,- Z;'^'''^'', with /jl'^'+'J being the nearest- 
neighbor two-body Hamiltonian density. Attached to each 
site is a three-indices tensor or Yg and to each bond 
a diagonal (singular value) matrix Aa or Ab, depending on the 
evenness and oddness of the /-th site and the /-th bond, respec- 
tively. Here, 5 is a physical index, s - 1 , ■ • ■ ,d, with d being 
the dimension of the local Hilbert space, and I and r denote 
the bond indices, l,r — \, - ■ ■ ,x, with x being the truncation 
dimension. The imaginary time evolution amounts to com- 
puting |»F(t)) = exp(-//T)|»F(0))/|exp(-//T)|»F(0))|. For large 
enough t and a generic initial state ^^(O)), it yields a good 
approximation to the ground state wave function, as long as 
there is a gap in the spectrum of the system. Following the 
Suzuki-Trotter decomposition the imaginary time evolu- 
tion operator is reduced to a product of two-site evolution op- 
erators acting on sites ; and U{i, - exp(-/!''''^'''5T), 
6t « 1. Notice that, a two-site gate U{i,i + I) renders the 
state not in the form of a MPS and breaks the translational 
invariance. The former is remedied by performing a singu- 
lar value decomposition of a matrix contracted from one F^ i^., 
one F^,^, one Aa and two Ab's, and only the x largest sin- 
gular values are retained. This yields the new tensors F^ 
F^ and Aa, which are used to update the tensors for all the 
sites, thus restoring the translational invariance under two site 
shifts. Repeating this procedure until the ground state energy 
converges, one may generate the systems's ground state wave 
functions in the MPS representations. 

Remarkably, for a system with symmetry-breaking or- 
ders, the iMPS algorithm automatically produces degenerate 
ground states arising from an SSB in the symmetry-broken 
phase, each of which breaks the symmetry of the system. 
Moreover, the symmetry breakdown results from the fact that 
an initial state has been chosen randomly. It is worth mention- 
ing that, for quantum lattice systems in one spatial dimension, 
continuous symmetries cannot be spontaneously broken Il23ll . 
due to strong quantum fluctuations [24]. Therefore, we shall 
restrict ourselves to the discussion of quantum lattice systems 
with a discrete symmetry group Z2 [25]. 

Bifurcations in the ground state fidelity per lattice site. Now 
consider a quantum many-body system, with a discrete sym- 
metry group Z2, on an infinite-size lattice in one spatial di- 
mension. Assume that the system undergoes a continuous 
QPT with Z2 symmetry spontaneously broken, when a con- 
trol parameter A varies. According to the definition JtI Htl, 
the ground state fidelity per lattice site, d{A\,A2), is the 
scaling parameter, which characterizes how fast the fidelity 
F{Ai,A2) = |(^(.l2)l*l'(^i))l between two ground states |>P(^i)> 
and '^{Ai)) goes to zero when the thermodynamic limit is ap- 
proached. In fact, the ground state fidelity F(Ai,A2) asymp- 
totically scales as F(Ai,A2) ~ d{Ai,A2)^, with L the num- 
ber of sites in a finite-size lattice. Remarkably, the ground 



state fidelity per lattice site is well defined in the thermody- 
namic limit, and satisfies the properties inherited from the fi- 
delity /I2): (i) normalization /I) = 1; (ii) symmetry 
d(AuA2) - diA2,Ai); and (iii) range < d{Ai,A2) < 1. 

In the Z2 symmetric phase, the ground state is non- 
degenerate, whereas in the Z2 symmetry-broken phase, two 
degenerate ground states arise. Now let us see what this im- 
plies for the ground state fidelity per lattice site, d(Ai,A2). If 
we choose *I'(/i2) as a reference state, with A2 in the Z2 sym- 
metric phase, then the ground state fidelity per lattice site, 
d(Ai,A2), cannot distinguish two degenerate ground states 
l*J'±('ii)) in the Z2 symmetry-broken phase. Here, |*I'+(/li)) = 
P\^-(Ai)), with P being the operation generating the symme- 
try group Z2. This follows from the fact that {^'(A2)\V+(Ai)} - 
mMPl'i'Mi)} = m^2)\'i'-Ui)), for any large but finite 
size L. However, if we choose *P(/l2) as a reference state, with 
A2 in the Z2 symmetry-broken phase, then d{Ai,A2) is able 
to distinguish two degenerate ground states. Therefore, for a 
given truncation dimension a pseudo phase transition point 
A^ manifests itself as a bifurcation point IJSil . An extrapola- 
tion Xo X - °° determines the critical point Ac- Therefore, the 
pinch point, first introduced in Refs. HSl as an intersection 
of two singular lines to characterize phase transition points, is 
identified as a bifurcation point. The significance of this re- 
sult lies in the fact that, on the one hand, this establishes the 
connection between the fidelity approach to QPTs and the sin- 
gularity theory; on the other hand, it is of practical importance 
since it makes possible to locate phase transition points with- 
out the need to compute the derivatives of the ground state 
fidelity per lattice site with respect to the control parameter 

In contrast, the von Neumann entropy, a bipartite entan- 
glement measure, fails to distinguish degenerate symmetry- 
breaking ground states. This is due to the fact that the von 
Neumann entropy is fully determined by the singular value 
matrices Aa and Ab, whereas all the information concerning 
an SSB is encoded in the tensors F^ and F^ 

Models. As an illustration, let us consider two quantum 
systems with the symmetry group Z2. The first is the quantum 
Ising model in a transverse magnetic field on an infinite-size 
lattice in one spatial dimension. It is described by the Hamil- 
tonian: 



(1) 



where S^^ (fl — x, z) are the Pauli spin operators of the /-th 
spin 1/2, and A is the transverse magnetic field. The model 
is invariant with respect to the operation: ^i'' — » -S^^^ for 
all the sites simultaneously, thus it enjoys the Z2 symmetry. 
As is well known, it undergoes a QPT, with a critical point at 
= 1 126]. 

The second is the spin 1/2 XYX model in an external mag- 
netic field, with the Hamiltonian: 

H=J] (^t'l^f 1 + A„5l'l5l'+'l + + hsfj , (2) 

where 5 J^^ (a — x, y, z) are the Pauli spin operators of the /-th 
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spin 1 12, Ay denotes the anisotropy in the internal spin space, 
and h is an external magnetic field. The model possesses a 
Z2 symmetry, generated by the operation: S^^^ -S]^^ and 
S^^ -S\!\ Note that A,. < 1 and A,, > 1 coiTespond 
to easy-plane and easy-axis behaviors, respectively. The or- 
dered phase in the easy-plane (easy-axis) case arises from an 
SSB along the x{y) direction, with a non-zero order parame- 
ter, i.e., the magnetization (5 t') ({S y^)) below the critical field 
he. Here we shall choose Ay = 0.25, for which it is critical at 
h - he, with he ~ 3.210(6) from the quantum Monte Carlo 



simulation 112711 . 

Simulation results. In Fig. [T] we present the probability 
mass function for the quantum Ising model in a transverse 
magnetic field in the Z2 symmetry-broken phase (A - 1/2). 
Suppose a random variable K follows the binomial distribu- 
tion with parameters n and p, then the probability of getting 
exactly k successes in n trials is given by the probability mass 
function: Pr (K ^ k) ^ C^^il - p)"-\ for A; = 0, 1, 2, ■ ■ ■ , «, 
where = n\/(k\{n - ky.) is the binomial coefficient. Here, 
by a success we mean that the order parameter (5''') is pos- 
itive. Our data are presented for both n = 20 and n = 40, 
with the truncation dimension to be 8. This confirms that 
the probability for getting the ground state with the positive 
order parameter (5 i'^) each simulation run is p - 1/2. The 
same pattern occurs for other choices of the truncation dimen- 
sion x- Actually this is true for any model with Z2 symmetry 
spontaneously broken. Therefore, our results demonstrate that 
an SSB occurs in classical simulations of quantum systems on 
an infinite-size lattice in the context of the iMPS algorithm. 
In contrast, algorithms that simulate finite-size lattice sys- 
tems are forbidden to produce degenerate symmetry-breaking 
ground states, since an SSB only occurs in the infinite-size 
(thermodynamic) limit. 

In Fig. |2]we plot the ground state fidelity per lattice site, 
d(Ai,A2), for the quantum Ising model in a transverse field. 
Here, the transverse magnetic field strength A is the control 
parameter. If we choose ^{A2) as a reference state, with A2 in 
the Z2 symmetry-broken phase (as shown here, A2 - 0.9), 
then d(Ai,A2) is able to distinguish two degenerate ground 
states, with a pseudo phase transition point A^ as a bifurcation 
point Il28il . The critical value /Ic = 1.00015 is determined from 
an extrapolation of the pseudo phase transition point A^ for the 
truncation dimension (see the inset in Fig.|2]i, which is quite 
close to the exact value 1. Therefore, the IMPS algorithm 
enables us to locate the transition point accurately from the 
computation of d(Ai,A2), with moderate computational cost. 
We stress that such a scaling for finite values of the trunca- 
tion dimension x has been discussed for the von Neumann 
entropy S. 

We have also presented d(hi,h2) for the spin 1/2 XYX 
model in an external magnetic field on an infinite-size lattice 
in Fig. [3] Here, the external magnetic field h is the control 
parameter. If we choose 'V{h2) as a reference state, with /12 in 
the Z2 symmetry -broken phase (as shown here, h2 - 3.2), then 
d{h\,h2) is able to distinguish two degenerate ground states, 
with a pseudo phase transition point h^ as a bifurcation point. 
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FIG. 1 : (color online) The probability mass function for the quan- 
tum Ising model in a transverse magnetic field in the Z2 symmetry- 
broken phase (A= \ /I). If a random variable K follows the binomial 
distribution with parameters n and p, then the probability of getting 
exactly k successes in n trials is given by the probability mass func- 
tion: Pr (K = k) = C*p*(l - py-\ for yt = 0, 1,2,-- - where 
C* = n\/(k\(n - k)\) is the binomial coefficient. Here, by a success 
we mean that the order parameter (S J.'') is positive. Our data are pre- 
sented for both n = 20 and n = 40, with ^ = 8. This confirms that 
the probability for getting the ground state with the positive order 
parameter (5'-') each simulation run is p = 1/2. The same pattern 
occurs for other choices of the truncation dimension x- Therefore, 
an SSB occurs in classical simulations of quantum systems on an 
infinite-size lattice in the context of the iMPS algorithm. 
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FIG. 2: (color online) Main: The ground state fidelity per lattice site, 
d(Ai, A2), for the quantum Ising model in a transverse field. Here, the 
transverse magnetic field strength A is the control parameter. If we 
choose *I'(/l2) as a reference state, with A2 in the Z2 symmetry-broken 
phase (as shown here, A2 = 0.9), then d{Ai,A2) is able to distinguish 
two degenerate ground states, with a pseudo phase transition point 
A^ as a bifurcation point. Inset: The critical point is determined from 
an extrapolation of the pseudo phase transition point /i^ for the trun- 
cation dimension^. The fitting function is A;^ = A^ + ax^'', where 
A, = 1.00015, with a = 0.31612 and b = 1.98565. This indicates 
that we are able to locate the transition point accurately, with moder- 
ate computational cost. The accuracy may be further improved if the 
truncation dimension x is increased. 
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FIG. 3: (color online) Main: The ground state fidelity per lattice site, 
d(h[,h2), for the spin 1/2 XYX model in an external magnetic field. 
Here, the external magnetic field h is the control parameter. If we 
choose T(/!2) as a reference state, with /i2 in the Z2 symmetry-broken 
phase (as shown here, /?2 = 3.2), then d(hi,h2) is able to distinguish 
two degenerate ground states, with a pseudo phase transition point 
as a bifurcation point. Inset: The critical point is determined 
from an extrapolation of the pseudo phase transition point for the 
truncation dimension The fitting function is = h^ + ax^'', where 
h, = 3.20471, with a = 0.02653 and b = 0.86003. 



The critical point he — 3.20471 is determined from an extrap- 
olation of the pseudo phase transition point /i^ for the trunca- 
tion dimension;^, as seen from the inset in Fig. [3] 

Summary. We have demonstrated that the IMPS algorithm 
may produce degenerate, symmetry-breaking, ground states 
arising from an SSB, each of which results from a randomly 
chosen initial state. It is shown that an SSB is reflected as 
a bifurcation in the ground state fidelity per lattice site for 
quantum lattice systems undergoing QPTs with symmetry- 
breaking orders. Conceptually, this establishes the connec- 
tion between the fidelity approach to QPTs and the singularity 
theory. Practically, it is also important since it makes possi- 
ble to locate transition points without the need to compute the 
derivatives of the ground state fidelity per lattice site with re- 
spect to the control parameter, which is usually a formidable 
task. We illustrated the general scheme in terms of the quan- 
tum Ising model in a transverse magnetic field and the spin 
1 /2 XYX model in an external magnetic field on an infinite- 
size lattice in one spatial dimension. 

Finally, we point out that we may extend our investigation 
to quantum lattice systems in two and higher spatial dimen- 
sions, which requires the iPEPS algorithm |18]. This is cur- 
rently under active investigation |19]. 
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